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Abstract 

We present preliminary measurements of direct and indirect CP- violating asymmetries in the decays 
of neutral B mesons to ir'^n'Tr^ and X^vr^vr'^ final states dominated by the resonance, using a 
time-dependent maximum likelihood analysis. The data sample comprises 88 x 10^ T(45) — > BB 
decays (80.8 fb~^) collected with the BABAR detector at the PEP-II asymmetric-energy B Factory 
at SLAC. For the CP violation parameters, we measure 

Agf^ = 0.19 ±0.14 (stat) ±0.11 (syst), Agp = -0.22 ±0.08 (stat) ±0.07 (syst), 
Cp^ = 0.45 toil (stat) ±0.09 (syst), Sp^ = 0.16 ±0.25 (stat) ±0.07 (syst). 

For the other parameters in the description of the B^(B^) pir decay-time dependence, we obtain 

ACp^ = 0.38 toil (stat) ±0.11 (syst), A5p^ = 0.15 ±0.26 (stat) ±0.05 (syst). 
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1 Introduction 



In the Standard Model, CP-violating effects arise from a single complex phase in the three- 
generation CKM quark- mixing matrix Q. One of the central questions in particle physics is 
whether this mechanism is sufficient to explain the pattern of CP violation observed in nature. 
Recent measurements of the parameter sin2/3 by the BABAR |2|, ^ and Belle [Q, ^] Collaborations 
establish that CP symmetry is violated in the neutral S-meson system. In addition, these two 
experiments have studied CP-violating asymmetries in B decays to the charmless two-body final 
states iC+7r~ and tt+tt" [P, ^. The time-dependent asymmetry in tt+tt" is related to the angle a 
of the unitarity triangle. 

In this paper, we investigate CP violation using charmless B^/B^ decays to ■K~^iT'''ir^ and 
K^TT^TT^ dominated by the p^h^ intermediate state, where h = tt ov K. As in the case of 7r"^7r~, 
the pTT mode provides a probe of both direct CP violation and CP violation in the interference 
between mixing and decay amplitudes. The latter type of CP violation is related to the angle a. 
In contrast to vr'^vr", p^vr^ is not a CP eigenstate and four configurations {B^{B^) p^n^) have 
to be considered. Although this leads to a more complex analysis it benefits from a higher 
branching fraction (20 - 30 x 10"^) |, |g. 

The p resonance is broad (150 MeV/c^) and the p^TT^ state may receive contributions at the 
amplitude level from other decay channels {e.g., B^ p'^ir^). For this analysis, we restrict 
ourselves to the two regions of the H^tt^tt^ Dalitz plot dominated by ph and assign a label, p~^h~ 
or p~h^ , to each event depending on the kinematics of the /i^vr^vr" final state. In the following, 
we will use the p~^h~ or p~h^ labels with the above meaning. 

Defining At = tp^ — ttag as the time interval between the decay of B^^ and that of the other B^ 
meson in the event, Sji'ag, the decay rate distributions can be written as [11| 



4r 



/So::,(Ao=(i±^&)^,^ 



|At|/T 



1 + [{Sph ± ^Sph) sin(AmdAt) - {Cph ± ACp,,) cos(AmdAt) 
1 - ({Sph ± ^Sph) sm{AmdAt) - {Cph ± ACph) cos(AmdAt) 



(1) 



The time-integrated charge asymmetries A(Jp and measure direct CP violation. The time 
dependence is described by four additional parameters. In the case of the self-tagging pK mode, the 
values of these four parameters are known to be CpK = 0, ACpK = — 1, SpK = 0, and ASpK = 0. 
For the pvr mode, they allow us to probe CP violation. Summing over the p charge in Eq. |l|, 
and neglecting the charge asymmetry A(Jp, one obtains the simplified CP asymmetry between the 
number of B'^ and B^ tags, given by 

A^o/BO = (Nbo - Npo)/{NBO + Npo) ~ Sp^ sin(AmrfAt) - Cp^ cos(AmdAt) . (2) 

The parameter Cpj^ describes the time-dependent direct CP violation and 5'^,^ measures CP violation 
in the interference between mixing and decay related to the angle a. 

The parameters AC^tt and ASpj^ are insensitive to CP violation. The asymmetry between 
iV(SO^ ^ p+TT-) + N{B%^ ^ p-7r+) and ^¥(5°^ ^ p-Tr+) + N{B^p^ ^ p+vr") is described by 
ACpTT) while AS'p^r is sensitive to the strong phase difference between the amplitudes contributing 
to B^ — > pvr decays. The naive factorization model |^ predicts ACp^ ~ 0.4 while there is no 
prediction for AS" 
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The measurements of the six parameters A'^, A'^p, Cpj^^ ^Cp.,^, Spj^, and ASpj^ reported 
here are performed using events collected by the BABAR detector at the PEP-II asymmetric- 
energy B Factory between January 2000 and June 2002. This sample corresponds to an integrated 
luminosity of 80.8 fb~^ taken at the T{4S) resonance ( "on-resonance" ) , which represents 88 million 
BB pairs, and 9.6 fb~^ taken around 40 MeV below the resonance ("off-resonance"). 

We extract the yields and CP parameters using a time-dependent maximum likelihood analysis 
based on Eq. |l[ This paper is organized as follows: the BABAR detector is described briefly in 
Sec. |2|. The event reconstruction and selection procedure is given in Sec. ^. i?-related backgrounds 
and their treatment in the likelihood analysis are described in Sec. |^. The full maximum likelihood 
fit is discussed in Sec. |5|. Finally, the results and the evaluation of systematic uncertainties are 
given in Sees. ^ and |^. 



2 The BaBar Detector 



A detailed description of the BABAR detector can be found in Ref. |12|. Charged particle momenta 
are measured in a tracking system consisting of a 5-layer double-sided silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) filled with a gas mixture based on helium and isobutane. The 
SVT and DCH operate within a 1.5-T superconducting solenoidal magnet. The typical decay vertex 
resolution is around 65 fim along the beam direction for the fully reconstructed B^j^, and around 100 
to 150 fim for the partially reconstructed tagging -B{^ag- Photons are detected in an electromagnetic 
calorimeter (EMC) consisting of 6580 CsI(Tl) crystals arranged in barrel and forward end-cap 
sub-detectors. The vr^ mass resolution is on average 7MeV/c^. The flux return for the solenoid is 
composed of multiple layers of iron and resistive plate chambers for the identification of muons and 
long-lived neutral hadrons. Tracks from the signal B decay are identified as pions or kaons by the 
Cherenkov angle 6ch measured with a detector of internally reflected Cherenkov light (DIRC). The 
typical separation between pions and kaons varies from 8a at 2 GeV/c to 2.5cj at 4 GeV/c, where 
a is the average 9ch resolution. Lower momentum kaons are identified with a combination of 9ch 
(for momenta down to 0.7 GeV/c) and measurements of ionization energy loss, dE/dx, in the DCH 
and SVT. 



3 The Event Selection and Reconstruction 

Signal Bph candidates are reconstructed from combinations of two charged tracks and a vr'^ candi- 
date. The charged tracks are required to be inconsistent with being an electron based on dE/dx 
measurements, shower shape criteria in the EMC, and the ratio of shower energy and track mo- 
mentum. The photons from the vr'' must have an energy greater than 50 MeV, and a lateral shower 
profile variable |T3[| between 0.01 and 0.6. The invariant mass 771(77) of the photons must satisfy 
0.11 < 771(77) < 0-16 GeV/c^. Similarly, to form a p candidate, the invariant mass m{'K^TT^) of the 
charged track and vr'^ must satisfy 0.4 < m(7r^-7r'^) < 1.3 GeV/c^. If both the (vr^vr'') and (7r~7r°) 
pairs satisfy this requirement, the B candidate is rejected, as the 'k~^'k~t^^ might result from inter- 
fering p's, and cannot be associated with a definite p charge. The track used for the p candidate 
must be inconsistent with being a kaon based on dE/dx and DIRC information. Finally, we require 
|cos0,r| > 0.25, where is the angle between the charged pion in the rest frame of the p and the 
p fiight direction in the rest frame of the B. We refer to the track h in ph as the bachelor track. 
To reject two-body S-background, the invariant mass of the two charged tracks, and the invariant 
mass of the bachelor track and the vr'^ must be less than 5.14 GeV/c^. 
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Two kinematic variables, used in the maximum likelihood fit, allow discrimination of signal 
B candidates from fake B candidates due to random combinations of tracks and candidates. 
The first variable is the beam-energy substituted mass defined as 



where s is the square of the center-of-mass (CM) energy, Ei and pi are the total energy and three- 
momentum of the e~^e~ state in the laboratory frame, and pb is the three- momentum of the B 
candidate in the same frame. Signal events populate the tues region around the B mass with a peak 
resolution of around 2.6MeV/c^. Candidates are required to satisfy 5.23 < niEs < 5.29 GeV/c^. 
The second variable, AE = E*^ — \/s/2, is the difference between the reconstructed energy of 
the Bph candidate in the CM frame and the beam energy. The /S.E distribution for signal events 
with a pion bachelor track (pvr) peaks around zero, while the distribution for pK signal events 
with the vr mass hypothesis assigned to the true kaon track, is shifted by —45 MeV on average 
(the exact shift depends on the momentum of the kaon). Backgrounds from other decay modes of 
the B peak at different /S.E depending on the number of charged and neutral particles in the decay: 
two-body decays, three-body decays, and four-body decays peak at positive, approximately 0, 
and negative Ai?, respectively. In order to reduce background from other modes, we require 
-0.12 < < 0.15 GeV. 

The time difference At is obtained from the measured distance between the z positions (along 
the beam direction) of the B^^y^ and B^^^ decay vertices, and the known boost of the e"'"e~ system. 
The vertex of the B'^^^ is reconstructed from all tracks in the event except those from the -B^/j, 
and an iterative procedure ||2[ is used to remove tracks with a large contribution to the vertex y^. 
An additional constraint is obtained from the three-momentum and vertex position of the B^^^ 
candidate, and the average e"'"e~ interaction point and boost. We require |At| < 20 ps and a{At) < 
2.5 ps, where a{At) is the error on At estimated on a per-event basis. 

Discrimination between pn and pK signal events is accomplished using the Cherenkov angle 
measurement from the DIRC. Therefore, only ph candidates with bachelor track inside the geomet- 
rical acceptance of the DIRC are considered. The number of photons in the DIRC associated with 
the bachelor track must be greater than 5. In addition, the Cherenkov angle of the bachelor 
track is required to be inconsistent with the proton hypothesis. Finally, we reject events where the 
bachelor track is inconsistent by more than 4cj with both the pion and kaon hypotheses. 

Continuum qq (where q = u,d, s, c) events represent the dominant background source for charm- 
less B decays. To enhance discrimination between signal and continuum background, we use a neu- 
ral network (NN) that combines four discriminating variables: two kinematic variables related to 
the reconstructed p^ candidate, the p mass and cosd.,^, and two event shape variables, Lq = J2iPi 
and L2 = ^iPi X |cos(0y Jp, where p* is the momentum of track i belonging to the rest of the 
event in the CM frame and 9^^ ■ is the angle between the momentum of track i and the B thrust 
axis Tb in the CM frame. Optimization and training of the NN is performed using off-resonance 
data contained in the signal region, to reduce residual correlations of the NN with the kinematic 
variables used in the maximum likelihood fit. In addition, the Monte Carlo signal training sample, 
generated with a GEANT4-based Monte Carlo simulation [jl^ ], only consists of correctly recon- 
structed signal events to increase the discrimination against cross-feed from other decay modes of 
the B. The distributions of the NN output for correctly reconstructed pir events, pn events with a 
misreconstructed vr*^, and continuum background are shown in Fig. A cut is applied on the NN 
output at 0.54 in order to reduce the number of continuum events entering the likelihood analysis. 
In addition, the NN output is used as a discriminating variable in the likelihood. Approximately 
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Figure 1: Distributions of the NN output for correctly reconstructed pir signal, misreconstructed 
pn signal and continuum background. 



37% of the events have more than one candidate passing this selection. In this case, we choose the 
candidate with the reconstructed tt^ invariant mass closest to the nominal vr*^ mass |T5[] . 

After all selection criteria have been applied, a total of 21172 events enter the likelihood fit. 
For 86% of the neutral B decays to p^n^ that pass the event selection, the charge of the p is 
unambiguously determined by the charge of the slowest track. If the lower-momentum track has 
a momentum below 2.4 GeV/c, its charge is assigned to the p. This approach does not rely on 
the reconstruction of the vr*^, and hence provides a more robust way to assign the charge to the p 
than criteria based on the reconstructed mass of p candidates. For the remaining events, the sign 
of the p is that of the vr^vr'' combination with invariant mass closest to the p mass [15]. With this 
procedure, only 5% of the events are assigned an incorrect charge. 

To determine the flavor of the B^^^ meson we use the same i?-tagging algorithm used in the 
BABAR sin2/3 analysis ||2|. The algorithm relies on the correlation between the flavor of the b quark 
and the charge of the remaining tracks in the event after removal of the tracks from the B ^ ph 
candidate. We define five mutually exclusive tagging categories: Lepton, Kaon, NTl, NT2, and 
Untagged. Lepton tags rely on primary electrons and muons from semileptonic B decays, while 
Kaon tags exploit the correlation in the process 6 ^ c — > s between the net kaon charge and 
the charge of the b quark. The NTl (more certain tags) and NT2 (less certain tags) categories are 
derived from a neural network that is sensitive to charge correlations between the parent B and 
unidentified leptons and kaons, soft pions, or the charge and momentum of the track with the 
highest CM momentum. The addition of Untagged events provides a larger sample for measuring 
the charge asymmetries and A^Jp. 

The quality of tagging is expressed in terms of the effective efficiency Q = J2c^c {Dc)'^, where 
Cc is the fraction of events tagged in category c and the dilution (Dc) = 1 — 2 (wc) is related 
to the average mistag fraction (wc)- The mistag fraction (wc), the efficiency Ec, and the mistag 
difference Awc = Wc — Wc, where Wc, Wc are the mistag probabilities for B^ and B'^, are measured 
for each tagging category c with a large data sample of fully reconstructed neutral B decays to 
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Table 1: Tagging efficiency €c, average dilution (Dc) = 1 — 2 (wc), dilution difference ADc = —2Awc, 
and effective tagging efficiency Qc for signal events in each tagging category. The values are 
measured with fully reconstructed neutral B decays. 



Category 


ec{%) 


(Dc) (%) 


AD, (%) 


Qc (%) 


Lepton 


10.7 ±0.2 


83.8 ± 1.6 


-0.4 ±2.4 


7.5 ±0.3 


Kaon 


34.8 ±0.3 


66.0 ± 1.2 


3.6 ± 1.6 


15.1 ±0.5 


NTl 


7.7 ±0.2 


58.4 ± 2.4 


-1.4 ±3.6 


2.6 ±0.2 


NT2 


14.1 ±0.2 


25.4 ±2.0 


7.2 ±3.0 


0.9 ±0.2 


Untagged 


32.7 ±0.4 








Total Q 








26.2 ±0.7 



D^*^ x+ (x+ = Tr^,p'^,ai ) and J/'ipK*'^ {K*^ A'+vr ) flavor eigenstates [§]. The tagging quality 
factor Q is found to be (26.2 ± 0.7)% (see Table |). We use these measurements of the tagging 
efficiencies and dilutions for pn and pK signal, and we float separate continuum background event 
yields for each category in the maximum likelihood fit. 

4 S-related Backgrounds 

We use a Monte Carlo simulation to study the potential cross-feed from other charm and charmless 
B decays starting from a list of more than 80 charmless decay modes to two-body, three-body and 
four-body final states, and an inclusive Monte Carlo simulation of i? — > charm decays. We estimate 
the number of events passing the event selection criteria using the selection efficiency from Monte 
Carlo and either measured branching ratios [^] or upper limits where available, or estimates based 
on related measured decay modes. We identify the 20 charmless modes that have more than one 
event entering the final sample. 

These modes are grouped into seven classes for which the discriminating variables have similar 
distributions. For each of the seven classes, a correction term is introduced in the likelihood, with 
a fixed number of events. Two additional classes for B~^ charm and B^ — > charm decays are also 
included in the S-background model, which is summarized in Table 

Like the selection efficiencies, the shapes of the distributions of the discriminating variables are 
obtained from Monte Carlo simulations. Figure § shows the AE-niEs planes for three main B- 
related backgrounds: B^ p\+ , B^ p'^p , and B — > charm. The charmless -B-background NN 
output and thes distributions are signal-like, and the AE variable discriminates between two-body 
{AE > 0), three-body (AE peaking around 0) and four-body (AE < 0) modes. The mES and 
AE distributions for B charm background have shapes similar to the continuum distributions. 

For charged i3-backgrounds, the At distribution is modeled as 

(At) = ^[l±A,{l-2wt)]e-\^'\/\ 
g^^o' (At) = ^[lTA,(l-2u;±)]e-l^*l/^ (4) 

tag.c 4tT 

where is the asymmetry between the number of p^h~ and p^h^ candidates for a given flavor tag 
and is extracted from Monte Carlo, and is the mistag fraction for tagging category c, measured 
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in data using a sample of fully reconstructed charged B decays to 0*^11^ for which we assume 
direct CP conservation. 

For neutral i?-backgrounds, the At distribution is parametrized as for signal, where A'^p = Sh = 
ASh = Ch = 0, and ACh is computed from Monte Carlo to take into account possible correlations 
between the reconstructed p charge and the flavor tag. We do not model CP violation for the 
i3-background in the nominal fit. The corresponding systematic uncertainties in our measurements 
are discussed in Sec. ^. 

Table 2: i?-background modes retained in the maximum likelihood fit, classified into nine categories. 
The number of expected events contributing to the pir and pK PDFs, scaled to 80.8 fb~^, and 
integrated over the full fit region, is reported in the 2nd and 3rd columns. The 4th and 5th 
columns give the assumed parameter values for the At distributions (see text). The last column 
gives the branching ratio (in units of 10"^) if measured, or the estimated range if not (these cases 
are indicated by the symbol *). 



Charg 


;ed Mode 


TV" 

cxp 


cxp 




Ak 


(Br 


± error) (10"^) 




■ i^*+(K+7r")pU 


0.41 


2.16 


1 


-1 




3 - 13* 


B+ ^ 


► K*^{K+TT )p+ 


0.06 


7.23 


1 


-1 




10 - 40* 


B+ ^ 




15.17 





0.22 






10 - 20* 


B+ ^ 


> ri'{p^-i)K+ 


0.19 


7.78 


1 


-1 




22.1 ± 2.1 


B+ ^ 


^ 7?'(p07)7r+ 


1.53 





-1 






1 - 5* 






23.22 





-1 






7 - 23* 


B+^. 




2.05 


21.20 


0.80 


-1 




8.4 ± 4.0 


B+ ^. 


■ p^Ti+ 


36.51 





-0.46 






9.7 ± 3.2 


B+^ 




8.24 





-0.76 






8.7 ± 1.3 


B+^ 


> K-^foiir+TT-) 


1.69 


15.46 


1 


-1 




11.7 ± 4.0 


B+^ 




0.17 


13.50 





-1 




11.6 ± 1.5 


B+^. 




5.09 





-1 







5.9 ± 1.4 


B+^. 


• charm 


195.0 


31.8 





-0.46 






Neutral Mode 


^ 'cxp 


^ ' cxp 


AC, 


ACk 


(Br 


± error) (10"**) 


B^^ 




3.38 





1 






8.7 ± 3.0 


B^^ 


K*+{K+TT^)p- 


0.72 


5.57 


1 


-1 




- 20* 


B^^ 


p+p- 


87.64 












40 - 100* 




pOpO 


1.03 












0-3* 


B^^ 


a^(/7r+)^- 


9.43 












28 - 48* 









6.62 




-1 




0-6* 


B^^ 


K*+{K+tt'^)tt- 


20.92 


12.97 


0.85 


-1 




8.7 ± 3.0 




K+TT- 


1.87 


2.17 


1 


-1 




18.5 ± 1.0 


B^^ 


charm 


121.7 


13.7 













5 The Maximum Likelihood Fit 
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Figure 2: AE-tties plane for three B-related backgrounds (Monte Carlo simulation): 
i?*' ^ p~^p~ and i? ^ charm. The lines indicate the cut on AE applied in the analysis; we require 
—0.12 < AE < 0.15 GeV. It removes most of the four-body i?-background modes (such as 
— > p^p^) and B charm modes. 



5.1 The Likelihood 

The yields and the other CP and non-CP observables are determined by minimizing the quantity 
— 21n£, where C is the total extended likelihood defined over all tagging categories c as 

5 Nc 

C = Y[e-^'llC,,,, (5) 

C=l 1 = 1 

with A'^^, the number of events expected in category c and Ci^c is the likelihood computed for event i. 

The sample is assumed to consist of signal, continuum background and S-background compo- 
nents where the bachelor track can be a pion or a kaon. The variables rriEs, AE and NN output 
discriminate signal from background, while the Cherenkov angle 6ch and, to a lesser extent, AE' 
constrain the relative amount of pir and pK. The variable At allows the measurement of the pa- 
rameters in the description of the B^{B^) pvr decay-time dependence and provides additional 
background rejection. The likelihood £j.c for event i in tagging category c is the sum of the 
probability density functions (PDF) over all components, weighted by the expected yields for each 
component, 

where 

• Nph is the number of signal events of type ph in the entire sample {h = tt, K). 

• ec is the fraction of signal events that are tagged in category c (given in Table ||). 

• ^qp/i is the number of continuum background events with bachelor track of type h that are 
tagged in category c. 
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• Vph = VP^{mEs) ■ VP^{AE) ■ VP^iNN) ■ VP^iOch) ■ VP^{^t) is the PDF for signal events. 
VP^{At) contains the measured physics quantities defined in Eq. |l| diluted by the effects of 
mistagging and the At resolution. 

• 'Pfc PDF for continuum background events with bachelor track of type h. 

• Cc'^ ^^'^ ^c'^ 3-^6 the i?-background contributions where the bachelor track is a true pion 
and a true kaon, respectively (see Sec. |^. 

Due to the relatively large number of low-energy photon candidates in reconstruction, the decay 
p^TT^ — > 7r"'"7r~7r° can be misreconstructed. These misreconstructed events have different shapes 
than correctly reconstructed signal for the distributions of the variables ttles, AE and NN output. 
Additionally, in some cases the assignment of the p charge may be wrong. These effects are taken 
into account by splitting the signal PDFs into three parts: true signal events that are correctly 
reconstructed, misreconstructed true signal events with right-sign charge, and misreconstructed 
true signal events with wrong-sign charge. The fractions of the three species of signal events are 
extracted from the Monte Carlo. 

5.2 The Probability Density Functions 

• ruEs 

The distribution for correctly reconstructed signal is parametrized using a Gaussian with a 
power law tail on the low side, where the mean is free to vary in the likelihood fit. The 
continuum background is parametrized using an ARGUS function with a floating shape 
parameter. 

• AE 

The distribution for correctly reconstructed signal is parametrized using the sum of two 
Gaussians, while the distributions for misreconstructed signal (both for right-sign and wrong- 
sign p charges) are modeled with simple Gaussians. The mean of the core Gaussian for the 
correctly reconstructed signal is floated in the maximum likelihood fit, in order to be less 
sensitive to the energy calibration for vr*^. Continuum background is modeled by a linear 
function. 

• NN output 

The NN output PDFs for correctly reconstructed and for misreconstructed signal events 
are determined with the Monte Carlo. A small discrepancy is observed between the NN 
output distributions from Monte Carlo and from a data control sample of fully reconstructed 
— > D~p~^ decays, and is propagated to the systematic error. The continuum PDF for the 
NN output is determined with the off-resonance data. 

For the above three variables, the PDFs for correctly reconstructed p^h^ decays and for 

misreconstructed B^ — > p^h^ decays are obtained from Monte Carlo after applying all selection 
cuts. The parameters for the continuum PDFs are either determined from off-resonance data, or 
left free to vary in the final fit to the on-resonance data sample. 

• At 

The resolution function for correctly reconstructed and misreconstructed signal events is a 
sum of three Gaussians, identical to the one described in Ref. |^, with parameters determined 
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from a fit to a large data sample of fully reconstructed neutral B decays to D^*^~x^ (x'^ = 
■K~^ , , a^) and J/iIjK*^{K*^ K^it~). The continuum At distribution is parametrized 
using a triple Gaussian with a common mean and three distinct widths that scale the At 
per-event error. The six parameters describing the shape of the the At continuum PDF are 
free to vary in the maximum likelihood fit. For each tagging category and bachelor hypothesis 
(vr or K), a parameter similar to Ah for i3-related backgrounds (see Sec. is introduced to 
take into account the correlations between the charge of the p candidate and the flavor tag. 
The values of these parameters are determined using on-peak data. 

• Particle identification 

The identification of the bachelor track as a pion or a kaon is accomplished with the Cherenkov 
angle measurement from the DIRC. We construct two Gaussian PDFs from the difference 
between measured and expected values of Och for the pion or kaon hypothesis, normalized by 
the resolution. The DIRC performance is parametrized using a data sample of D*^ D^-k^ , 
K~tt'^ decays. Within the statistical precision of the control sample, we find a similar 
response for positive and negative tracks and use a single parametrization for both. 



6 Results 

We find 413I33 (stat) pir and 147121 (stat) pK events in our data sample and we measure the 
CP parameters 

^gf^ = 0.19 ± 0.14 (stat), A^Jp = -0.22 ± 0.08 (stat), 

Cp^ = OA5toil (stat), Sp^ = 0.16 ± 0.25 (stat). (7) 

The two other observables in the decay rates (Eq. |l|) are measured to be 

ACp^ = O.SSto.ll (stat), AS'p^ = 0.15 ±0.26 (stat). (8) 

The correlations between these parameters are summarized in Table ^. 

Alternatively, the results on direct CP violation can be expressed using the asymmetries 

, ^ N{B% ^ p+vr-) - iV(i?o, ^ p-^+) _ A"^p - Cp^ - A'^p ■ ACp^ 

+- N{BO^ - p+TT-) + N{B% - p-^+) 1 - ACp. - A^^p -Cp^ ' ^' 
, _ NjB^p^ ^ p-7r+) - NjBl ^ p+TT-) _ A^Jp + Cp^ + APJp.ACp^ 



N{B% - p-vr+) + N{B% - p+^-) 1 + ACp. + A'^p • C, 



p-TT 



In the decays i?^^ — > p^ir^ and Bp^ p~7r~^ the spectator quark is involved in the formation of 

the p meson. These two decay modes are related to the direct CP asymmetry A^ according to 

Eq. ^. Similarly in Eq. |l^, we probe direct CP violation through the asymmetry A |_ using the 

decays B^^ —>■ p~Tr~^ and Bp^ — > p+Tr". In this case the vr meson is formed from the spectator 
quark. From the above fitted values (Eq. |^ and Eq. |8|) and their correlation matrix (Table ^) we 
obtain 

= -0.82 ±0.31 (stat), = -0.11 ±0.16 (stat). 

Figure |^ shows the distributions of tues and AE for data samples that are enhanced in signal 
using cuts on the signal-to-continuum likelihood ratio of the other discriminating variables. Fig- 
ures |^ and |5| show respectively the time-dependent asymmetry Apo jpo (see Eq. ^) between B^^^^ and 
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Table 3: Correlation coefficients (in percent) between the six parameters A^p, Cpj^, AC/^, 

SpTT and ASpTT obtained for a data sample of 88 million BB pairs. The global correlation coefficient 
is the largest correlation between the parameter in question and any linear combination of the other 
28 free parameters in the likelihood fit. 





Global 
correlation 


^CP 


^CP 




ACf^ 




ASp^ 


^CP 


15.6 


100 


3.4 


-11.8 


-10.4 


0.6 


0.5 


^CP 


6.3 


3.4 


100 


-1.3 


-1.1 


-0.4 


-0.5 


C 


28.8 


-11.8 


-1.3 


100 


23.9 


9.2 


-6.8 


AC 


28.3 


-10.4 


-1.1 


23.9 


100 


6.9 


-9.2 


S 


24.9 


0.6 


-0.4 


9.2 


6.9 


100 


-23.4 


AS 


25.1 


0.5 


-0.5 


-6.8 


-9.2 


-23.4 


100 



-^tag events in the Lepton and Kaon categories, and the time-dependent asymmetry Ap+^- ip-^+ 
between p^tt~ and p~tt'^ for all the tagging categories, after a cut on the signal-to-continuum 
likelihood ratio of all discriminating variables except At. 

As a validation of the At parametrization in data, we allow r and Arrid to vary in the fit. We 
find r = (1.59 it 0.12) ps and Am^ = (0.51 it 0.09) ps^^, and the remaining free parameters are 
stable with respect to the nominal fit with fixed r and Am^- When allowing ACpK to vary, we find 
ACpK = —1.15 lb 0.34, which is consistent with the expected value of —1 used in the nominal fit. 
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Figure 3: Distributions of uies and AE for samples enhanced in p-zr signal using cuts on likelihood 
ratios. The solid curve represents a projection of the maximum likelihood fit result. The dashed 
curve represents the contribution from continuum events {pir and pK candidates combined) , and the 
dotted line indicates the combined contributions from continuum events and B-related backgrounds, 
including pK. 
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Figure 4: Time distribution and asymmetry (^^o /^o in Eq- @) for B^^^ and B'^^^ decaying to 
pvr, in the Lepton and Kaon categories. The sample was enriched in signal events by applying a 
cut on the signal-to-continuum likelihood ratio. The solid curve is a likelihood projection of the 
result of the full fit, and is normalized to the expected number of events according to that fit (71 
signal events, 36 continuum background events and 10 B background events). The dotted line is 
the contribution from i?-related backgrounds and the dashed line is the total B and continuum 
background contribution. The depression around zero in the asymmetry plot is due to the relatively 
large fraction of continuum events in this region of At. The non-zero central value for the CP 
parameter Spj^ induces a small contribution to the asymmetry that is odd in At. 
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Figure 5: Time distribution and asymmetry Ap+^- /p-^+ between p'^ir" and p^ir^ for all the tagging 
categories. The sample was enriched in signal events by applying a cut on the signal-to-continuum 
likelihood ratio. The solid curve is a likelihood projection of the result of the full fit, and is 
normalized to the expected number of events according to that fit (156 signal events, 157 continuum 
background events and 21 B background events). The dotted line is the contribution from B- 
related backgrounds and the dashed line is the total B and continuum background contribution. 
The depression around zero in the asymmetry plot is due to continuum dilution. In the absence of 
continuum background, the asymmetry curve would be flat and equal to A^Jp defined in Eq. |l[ 
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7 Systematic Uncertainties 



Several sources contribute to the systematic error on the measured observables. 

• A?n,(i and B lifetime 

This systematic error is estimated by varying the values of Am^ and the B lifetime with 



respect to their world average |15] by one standard deviation 



• At resolution function 

Since the parameters of the time distribution for the continuum background are fitted in 
the data, the only contribution to the systematic uncertainty is due to the At resolution for 
signal. We propagate the errors on each parameter (scale factors, biases, fraction of tails) of 
the signal At resolution function to the fit results. 

• Tagging 

We propagate the uncertainties on the tagging dilutions and the tagging efficiencies to the fit 
results. 

• Fraction of misreconstructed signal events and misreconstruction of p charge 

Both the the fraction of misreconstructed signal events and the fraction of events with mis- 
reconstructed p charge and are determined with the Monte Carlo. We assign a conservative 
25% uncertainty to these parameters and propagate this error to the fit results. 

• Fitting procedure 

We perform fits on large Monte Carlo samples of a weighted mixture of pn/pK signal, as well 
as fits to a large number of toy Monte Carlo samples of the same size as our data sample that 
contain realistic amounts of both continuum and i?-related backgrounds. The small biases 
observed in these two tests are added in quadrature and assigned as a systematic uncertainty. 

• ruES, and NN output PDFs 

The continuum PDFs for uies and AE are fitted to the data in the likelihood analysis. 
The continuum PDF for the NN output is determined with the off-resonance data and we 
propagate the statistical uncertainty due to the limited size of this sample to the fit results. 
We evaluate the systematic uncertainties due to the signal PDFs with a large B^ D~p^ 
data control sample. The small differences observed between the distribution shapes for 
Monte Carlo events and the distribution shapes obtained from the data control sample are 
used as an estimate of the systematic uncertainty on the signal PDFs. These uncertainties 
are propagated to the fit results. 

• Particle identification 

We assume that the average position of the Cherenkov angle {9ch) is known with a precision of 
0.5 mrad and that there is an uncertainty of 10% on cr{6ch)- We propagate these conservative 
errors to our fit results. 

• S-related backgrounds 

The number of events in the various modes entering our description of the B background are 
varied according to the error on their branching ratio, if measured, or in the range indicated 
in Table |. 

The parameters (A^ and ACh) describing the correlation between the tagging and the p charge 



21 



assignment are varied within a conservative range. For some of the major contaminations 
{p^p^, p^TT'^, and B — > charm) conservative ranges are extracted from various Monte Carlo 
studies. For other major contaminations (K*^ {K^ 71^)11^ , p^K) a conservative error equal 
to the correction itself is used. For the other modes, the full range [—1, +1] is used for the 
systematic study. 

All i?-background modes can potentially exhibit direct CP violation, and a few of them can 
potentially exhibit CP violation in the interference between decay and mixing. Finally, for the 
neutral modes, various physical phases may lead to a non-zero value of the phenomenological 
parameter A5/i. For the major contaminations listed previously, a Monte Carlo study yields 
the maximum possible range for the corresponding effective parameters. This range is then 
used for the systematic study. For the other modes, the full range [— is used. 

Table ^ summarizes the various sources contributing to the systematic error on the measure- 
ments of the six parameters ^ct', CpT^, ACpj^, Sp.,^, and ASp,r- The main source of systematic 
error arises from the uncertainty on the B background components. Finally, the systematic errors 

for the direct CP asymmetries Aj^ and A |_ (see definitions in Eq. |^ and Eq. [lO|) are ±0.16 and 

±0.09 respectively. 



Table 4: Summary of the systematic uncertainties on the CP and non-CP observables. The indi- 
vidual systematic errors are added in quadrature. 



Type of systematic error 


^CP 


^CP 


Cp-ji 


ACpTT 






Am ± 0.008 ps"^ 


0.000 


0.000 


0.005 


0.005 


0.001 


0.001 


r± 0.016 ps 


0.001 


0.000 


0.002 


0.002 


0.002 


0.002 


Time resolution 


0.001 


0.001 


0.003 


0.002 


0.005 


0.006 


Tagging 


0.001 


0.001 


0.021 


0.012 


0.012 


0.012 


Fraction of misrec. signal events (±25%) 


0.002 


0.002 


0.017 


0.004 


0.002 


0.004 


Fraction of events with misrec. p charge (±25%) 


0.000 


0.002 


0.001 


0.004 


0.000 


0.000 


Fitting procedure (bias) 


0.012 


0.001 


0.038 


0.033 


0.020 


0.015 


NN output signal {B^ ^ p+) 


0.002 


0.002 


0.006 


0.006 


0.009 


0.013 


NN output continuum (off-resonance data) 


0.002 


0.002 


0.001 


0.001 


0.001 


0.001 


PDF for ruEs and 


0.006 


0.003 


0.009 


0.003 


0.005 


0.017 


Particle identification 


0.010 


0.010 


0.008 


0.014 


0.013 


0.001 


^-backgrounds 


0.109 


0.065 


0.077 


0.107 


0.063 


0.035 


Total 


0.110 


0.065 


0.091 


0.114 


0.069 


0.046 



8 Summary 

With a data sample of 88 million BB pairs, collected between January 2000 and June 2002 by the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC, we find 413^13 (stat) p-rr and 
147^21 (stat) pK events and we obtain the following preliminary measurements of the CP violation 
parameters: 

= 0.19 ± 0.14 (stat) ±0.11 (syst), A^Jp = -0.22 ±0.08 (stat) ±0.07 (syst), 
Cp^ = 0.45 toil (stat) ±0.09 (syst), Sp^ = 0.16 ±0.25 (stat) ±0.07 (syst), 
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and of the other parameters in the description of the B^{B^) pir decay-time dependence: 
ACp^ = 0.38 toil (stat) ±0.11 (syst), ASp^ = 0.15 ±0.26 (stat) ±0.05 (syst). 

For the asymmetries j4+_ and A |_, which probe direct CP violation, we measure 

= -0.82 ±0.31 (stat) ±0.16 (syst), A_+ = -0.11 ±0.16 (stat) ±0.09 (syst). 
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